Meteor shower identification and characterization 
with Python 

Althea Moorhead 
Meteoroid Environment Office 
NASA Marshall Space Flight Center- 






All Sky Fireball Network 
Shower Identification with Python 
Other MEO Python Applications 


Outline 


Q All Sky Fireball Network 
o Shower Identification with Python 
o Other MEO Python Applications 


NASA’s Meteoroid Environment Office (MEO) is the 
NASA organization responsible for meteoroid 
environments pertaining to spacecraft engineering and 
operations. 
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Fireball Networks 



o NASA All Sky Fireball Network 
o Southern Ontario Meteor Network (SOMN - UWO) 
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Events 



o Overlapping fields of view for trajectory triangulation 

© Automatic meteor detection using ASGARD 
(Weryk et al. 2008) 


Althea Moorhead, NASA MEO 


Meteor showers with Python 


4/ 15 




All Sky Fireball Network 
Shower Identification with Python 
Other MEO Python Applications 



o Radiant, speed, and atmospheric height 
o Geocentric trajectory and heliocentric orbit 
© Categorized using shower surveys (IAU MDC, Brown et 
al. 2008 ) 
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Shower identification using orbit 
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angles 


9 Need agnostic approach 
for new/outbursting 
showers 

9 Need thorough shower 
removal 

o Showers show up clearly 
in orbit angles, 
sometimes more 
obviously than in 
radiant space (KCGs, 
ZCSs) 
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Shower membership using orbital similarity 


o Drummond (1981) D parameter: 



o Is often used for both parent body identification and meteor 
classification 

o We Monte-Carlo meteor orbits using uncertainties and look 
for overlap with showers 
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Algorithm 


9 Python code: 

a. Bin meteors by date, /, and ft 

b Take most populated bin and compute average orbit. 

a Weight meteors by 1/D and iterate, 
c Compare to list of showers. If there is no match, output bin 
meteors and exit. 

d Assess shower membership using D for each meteor, taking 
uncertainties into account (5% overlap counts as membership) 
e. Remove shower meteors, go to a. 

Q Inspect for similarity to known shower 

Q Add shower to list of showers. Go to 1. 
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2014 k , Cygnid outburst 
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Work in progress: D weighting 


o Detects major showers 
(and outbursts!) 

o Identified minor showers that 
made significant contributions 
(e.g., ZCS) 

9 Ongoing: 

o Use D to detect orbit 
clusters 

o Quantify shower 

membership probability 
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Atmospheric trajectories of meteors 
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Atmospheric trajectories of meteors 


Uses Scipy’s ODE modules: 


# Calculate trajectory through atmosphere from orbit 

def get_traj (xv0, pm f cd, dt=0.1, ablate=False, frag=False) : 

# Do the integration, stop with landfall or escape 

r = ode(derivs) . set_integrator( 1 v ode 1 , method=' bdf ' , 
with_jacobian=False) 
r.set_initial_value(xv0, 0.0) 
r.set_f_params(pm, cd, ablate) 

# Integrate while in atmosphere and at least a cm in diameter 
ts, Is, ps, hs, xvs , dt = [], [], [], [], [], 0.1 

sz = r.y[0] 

while r. successful( ) and (0 < h_r(r.y) < l.l*h_atmo) and (r.y[0] > minsz) : 

if r.y[0] - sz > 1.0e-10 : 
print sz, r.y[0] 

raise Exception ( 'meteoroid has grown!') 
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Annual shower forecasts 
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Annual shower forecasts 


2015 Fluxes 

— 0.04 cm 0.30 cm 



UT Date 


o Translated from IDL to Python, using array broadcasting to 
eliminate FOR loops 

o New version is significantly faster 
o Python version of VSOP87 solar longitude calculation 
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Comet/meteoroid ejection model 


Q JPLEphem for comet 
position 

Q Gravity + radiation 
pressure: 

P GMm „ 

F = -{\-p) r 

O Solve Kepler’s equation 
(Gooding & Odell, 1988) 

— 10 6 particles in seconds. 
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